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ABSTRACT 

We present soft X-ray variability time-scales for 65 active galactic nuclei 
(AGNs) derived from ROSAT/PSPC pointing data. For these 65 objects with 
obvious exponential increasing or decreasing patterns in their light curves, we use 
the exponential time-scales and find they are more suitable for analyzing ROSAT 
light curve data. We also discuss some physical constraints on the central engine 
implied by our results. Assuming that this soft X-ray variability exponential 
time-scale is approximately equal to the thermal timescale of the standard ac- 
cretion disk, we obtain the accretion rate, the size of soft x-ray radiation region 
and the compactness parameter for 37 AGNs, using their recently estimated cen- 
tral black hole masses. For 12 of these 37 AGNs, the radii of the gravitational 
instability in the standard thin accretion disks are obtained using central black 
hole masses and the calculated accretion rates. These are consistent with the re- 
sults from the reverberation mapping method. These results provide supporting 
evidence that such gravitational instabilities contribute to the formation of the 
Broad Line Regions (BLRs) in AGNs. 

Subject headings: X-rays: galaxies — accretion, accretion disks — galaxies: ac- 
tive — galaxies: nuclei — galaxies: Seyfert — quasars: general 



1. INTRODUCTION 

X-ray variability has long been known to be a common property of active galactic nuclei 
(AGNs). The X-ray variability in the continuum and in the emission lines can give a clue to 
the central engine in AGNs (Mushotzky et al. 1993). The variance amplitude and the time- 
scale are two main parameters which characterize the variance. There are many quantities 
can be used to describe the variability amplitude or the timescale in X-ray variability, such 
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as the flux-doubling time-scale (Barr & Mushotzky 1986), normalized variability amplitude 
(Green et al. 1993), and the excess variance (Nandra et al. 1997). Some powerful methods 
are also introduced in the analysis of the X-ray light curves. These include the power density 
spectrum (PDS), the wavelet method (Otazu et al. 2002 ), and the power density in the 
time domain (Li 2001). 

There has been steady progress on the analysis of X-ray variability since the launch 
of EXOSAT, RXTE, ROSAT, ASCA, XMM-Newton. For an individual X-ray light curve, 
PDS analysis is now popular and can offer more information about the X-ray variability. 
Early attempts to constrain a PDS flattening or "cut off", which relates to a characteristic 
time-scale, were made with the EXOSAT data of NGC 5506 (McHardy 1988) and NGC 
4151 (Papadakis 1995). These analyses yield evidence for a cutoff time-scale of several 
weeks. However, the uneven sampling of these data made their reliability uncertain. The 
situation was significantly improved with the launch of RXTE. Using RXTE data, Edelson 
& Nandra (1999) obtained an evenly sampled X-ray light curve of NGC 3516 and produced 
a PDS covering four decades in temporal frequency, finding a cutoff time-scale of about one 
month. Pounds et al. (2001) also found a cutoff time-scale of about 13 days for a narrow 
line Seyfert 1 galaxy (Akn 564). Uttley et al. (2002) recently developed a Monte Carlo 
method to test models for the true power spectrum of intermittently sampled light curves 
against the observed noisy power spectra of four Seyfert galaxies, which are flat towards lower 
frequencies. They fit two models for the flattening: a "knee" model and a "high-frequency 
break" model. They reduced the characteristic frequency for these four AGNs. 

The X-ray characteristic time-scale or frequency usually relates to the size of the X-ray 
emission region, which then can be used to constrain the central engine, e.g. the central 
black hole masses, accretion rates, and the radiation mechanism. 

ROSAT observations of AGNs offer plenty of information about the soft X-ray vari- 
ability, which is thought to come from the accretion disk around the central black hole in 
these AGNs. However, PDS analysis is not suitable for analyzing ROSAT pointed data, 
since there are usually large time gaps between adjacent exposures in ROSAT pointing ob- 
servations. Thus, time-scale analysis remains necessary in AGN samples. There are obvious 
exponentially increasing or decreasing patterns in the soft X-ray light curves in some AGNs 
observed with ROSAT/PSPC in pointing mode. In this paper we use the exponential time- 
scale (Zhao & Fink 1996) to investigate the soft X-ray variability in AGNs with this kind of 
pattern. This paper is a continuation of the previous work of Zhao & Fink (1996) . 

The standard model of AGNs is an accretion disk with a jet around a central massive 
black hole. Progress on a reliable estimate about black hole masses can effectively constrain 
the central engine when combined with the research on the X-ray variability. Several methods 
have been used to estimate the mass of a black hole in the center of a galaxy. These include 
stellar dynamical studies (reviewed by Kormendy & Gebhardt 2001) , reverberation mapping 
methods (Kaspi et al. 2000), analyses of the high frequency tail of the power density spectrum 
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(Czerny et al. 2001), the relation between central black hole masses and the bulge velocity 
dispersion Mbh—o (Merritt & Ferrarese 2001), and the use of single-epoch rest-frame optical 
spectrophotometric measurements (Vestergaard 2002). 

In next section time-scales and central black hole masses of our sample are presented. We 
then give the physical constraints on the central engine of 38 AGNs in Section 3. In Section 
4 we derive the radii of the gravitational instability in the standard a accretion disks of 12 
AGNs and compare them with BLRs sizes from the reverberation mapping method. Our 
results are summarized in Section 5. The cosmological parameters H = 75kms _1 Mpc _1 and 
go = 0.5 have been adopted in this work. 



2. SAMPLE AND TIMESCALES OF VARIABILITY 

All sources were observed in ROSAT/PSPC mode over periods varying from days to 
years. We picked out the AGNs from cross identification of the Veron-Cetty (2001) AGN 
catalogue with the ROSAT point source catalog. From the ROSAT public archive of PSPC 
observations, only sources with total X-ray photon counts greater than 1000 were selected 
to ensure the quality of the X-ray spectra. This yielded more than 200 AGNs. The data 
were processed for instrument corrections (such as vignetting and dead time effects) and 
background subtraction using the EXS AS/MIDAS software. 

The light curve for each AGN was obtained from original ROSAT observations with a 
400 seconds time bin in three energy bands; 0.1-2.4 (total band), 0.1-0.4 (A band), 0.5-2.0 (B 
band) keV. The estimates of time-scales for all AGNs described below are based on those light 
curves. From these light curves we can find that there exist obvious exponentially increasing 
or decreasing patterns. Since X-ray emission is believed to come from the accretion disk 
around the central black hole in AGNs, these exponential patterns may provide physical 
insight. We select the AGN light curves with exponential patterns by the following criteria: 
(l)The count rate increased or decayed continually in the light curves, considering count rate 
uncertainties of 3 a, (2)the number of the contiguous data points fitted by the exponential 
function (Eq. 1) are greater than 4, (3)the magnitude of the count rate variance is greater 
than 5a/ (of + of), where o\ and 02 are the count rate errors of the jumping-off point and 
end-point, respectively. The soft X-ray luminosity in the 0.1-2.4 keV of each AGN is obtained 
by averaging the power-law fitting spectra of all the observations. 



2.1. Exponential Time-scales 

We use the exponential time-scale (Zhao & Fink 1996) to characterize the exponential 
pattern in the light curves in AGNs. The exponential time-scale, At e , is defined as 

I(t)=I + I a e ± ?&. (1) 
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where I(t) is the count rate at t, Iq and I a are constants. 

The definition of exponential time-scale is based on instabilities which increase or de- 
crease exponentially. Thus the exponential time-scale is equivalent to a rising or a descending 
time-scale and is the amplitude of variability. Iq is an additional constant component, 
which is equal to the intensity before or after the instability. It is exponentially increasing if 
the sign of I a is positive while it is exponentially decaying if the sign is negative. For AGNs 
with some values of At e , we select At e with the smallest error as the time-scale. Taking 
the Seyfert 1 galaxy NGC 5548 as an example (Fig. 1), there are three obvious exponen- 
tially increasing or decreasing patterns in the light curve. Three patterns fitting by Eq. 1 
showed similar exponential time-scales. The exponential time-scales are (1.04 ± 0.32) x 10 5 , 
(6.65 ±5.35) x 10 4 , (8.17±0.26) x 10 4 seconds, respectively, which are consistent if taking the 
errors into account. And I are 3.10 ±0.13, 7. 16 ±0.87, 6.01 ±0.10 photons s -1 , respectively. 

We think there is more physical information in the exponential time-scale than in the 
two-folding time-scale, which only considers the amplitude increased by a factor of two. 
There are usually large time gaps between adjacent exposures in ROSAT pointed observa- 
tions. The two-folding time-scale will overestimate the time-scale of ROSAT sources. Using 
our selection criteria as defined above, we searched for any obvious continually exponentially 
increasing or decreasing pattern in the light curve of each AGN and then fitted the pattern 
with the form of Eq. 1, considering the errors of the count rates (Press et al. 1992). The 
exponential time-scales fir 65 AGNs are listed in column 3 of Table 1, in which there are 
nine QSOs, 37 Seyfert 1 galaxies, nine NLSls, one Seyfert 2 galaxies, four BL Lacs and five 
high optical polarization (HP) objects. The time-scales calculated as rising exponentials are 
marked with a star. The rest are derived from decaying exponentials. 



2.2. Estimations of Central Black Hole Masses 

The masses of 13 AGNs labelled with "a" in Table 1 are adopted from Nelson (2000). 
They were estimated from the reverberation mapping methods with a high level of reliability. 
There is strong correlation between the black hole mass and the host velocity dispersion, 
namely the M BH — a relation (Ferrarese et al. 2001). We adopt the M BH — a relation found 
by Merritt & Ferrarese (2001), namely, 

M BH = 1.3 x lO 8 M (o-/2OOkms~ 1 ) 4 - 72 . (2) 

We derive masses for seven AGNs with measured host velocity dispersions, which are 
labelled with "b" in Table 1 (Nelson & Whittle 1995; Falomo et al. 2002). The errors for 
black hole masses in these seven AGNs are estimated from the error of the bulge velocity 
dispersion. An empirical relationship between single-epoch rest-frame optical spectropho- 
tometric measurements and the central masses was recently presented (Vestergaard 2002). 
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These methods for estimating AGN masses are all calibrated using masses from the reverber- 
ation mapping method. The total number of AGNs with available black hole masses is 37, of 
which there are 22 Seyfert 1 galaxies, five NLSls, one Seyfert 2 galaxies, eight Quasars and 
one BL Lac objects. The central black hole masses of these 37 AGNs are listed in column 5 
in Table 1. 



3. ANALYSIS AND PHYSICAL CONSTRAINTS 



3.1. Central Black Hole Mass and Mass Limits from the exponential 

Time-scale 

The upper limit on the mass of the central black hole in AGNs derived from the time- 
scale is defined as (Mushotzky et al. 1993): 

M(At) = 2 x 1O 4 (^-)M . (3) 

assuming that the size of the X-ray radiation region is 5 Schwarzschild radii. 

We compared this mass limit derived from the exponential time-scales with the above 
estimated central black hole masses (see Sec. 2.2) in Figure 2. This shows that the expo- 
nential time-scale can also provide an effective limit for the mass of a central black hole in 
an AGN. 



3.2. Constrains on the Radiation Mechanism in AGNs 

Some physical constraints on the relationship of bolometric luminosity and the time- 
scale were proposed by Fabian (1979; 1992). If the source luminosity is below the Eddington 
limit, we obtained a limit: 

L b oi < -n^—At ergs s \ (4) 

6 <Jt 

When the luminosity is produced by matter being transferred into radiation with some 
efficiency 77 <§C 1, there is a second limit: 

4 

L bo i < r]—?—At ergs s~\ (5) 

If the primary spectrum of the source extends to 7-ray energies, electron-positron pair pro- 
duction gives a third limit: 

4 

Lboi < 20tt^— At ergs s~\ (6) 
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Thermal bremsstrahlung radiation supplies more power than Compton process only if: 

L bol < 10 39 T 9 3/2 At ergs s" 1 . (7) 

where the temperature of the electrons is T 9 10 9 K. 

The bolometric luminosity is larger than the X-ray luminosity in ROSAT's energy band. 
For example, the ratio of the bolometric luminosity of 3C 273 (Padovani & Rafanelli 1988), 
1.7 x 10 47 ergs s _1 , and the X-ray luminosity in ROSAT's energy band, 1.2 x 10 46 ergs s _1 , 
is 14. In Figure 3 we show the soft X-ray luminosity in 0.1-2.4 keV versus the exponential 
time-scale. The lines in Figure 3 are from Eq. 4-7 assuming L bol = 10L X . 

From Figure 3 bremsstrahlung is ruled out as the primary radiation source except in 
several low luminosity sources because of the observed rapid variability. 

The objects whose bolometric luminosity are larger than the Eddington limit are MS 
01585+0019 (BL) ,PKS 0548-322 (HP) ,MS 23409-1511 (NLS1), PG 1244+026 (NLS1) . 
These objects may display relativistic effects. 

The efficiency factors, 77, of most objects are larger than 0.007, which would not be 
associated with stellar process (see Fig. 3) . It suggests strongly that most AGNs, at least, 
have non-thermal and non-stellar processes which are produced by a very efficient central 
engine. 

We also plot L x versus central black hole masses in Figure 4. The solid line shows 
L x = O.lLEdd- and the dash line is L x = 10~ 4 LEdd- If we assume L bo \ = 10L X , we find that 
the efficiency factor is between 1 to 10 -3 , which is consistent with the results from L x — At 
diagram ( see Fig. 3). 



3.3. Sizes of the X-ray Radiation Regions 

Under unified model schemes, an accretion disk surrounds the central massive black hole 
in an AGN. There are two instabilities, thermal and viscous, usually occurring in accretion 
disks. However, the rapid variability in X-rays cannot arise from viscous instability because 
the viscous time-scale is too long (Mushotzky et al. 1993). The time-scale of thermal 
instability in accretion disks (Frank et al. 1992) is 



where a is a parameter in the standard thin accretion disk, Q k = y ^+Jf is the Keplerian 
angular velocity in the orbit of radius R surrounding a central black hole with the mass of 
M. 
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Assuming the exponential time-scale is approximately equal to the thermal time-scale, 
we can obtained the size of radiation region, R x , as 

Rg \pl Rg 

where R g = is the Schwarzschild radius, a — 1 is adopted. 

Taking the above estimated black hole masses and exponential time-scales, we can ob- 
tained the ratios of the size of X-ray radiation region to the Schwarzschild radius, which are 
listed in column 7 of Table 1. R x /R g distribution is < log(R x / R g ) >= 0.97±0. 12 with a stan- 
dard deviation of 0.75. We plot R x /R g versus Mbh in Figure 5. A simple least square linear 
regression (Press et al. 1992) gives log(R x /R g ) = 4.52±0.92 + (-0A7±0.12)log(M BH / M ), 
with a Pearson correlation coefficient of R=0.54 corresponding to a probability of P<0. 000431 
that the correlation is caused by a random factor. The error of R x /R g is estimated from the 
errors of the exponential time-scale and the black hole mass. The sizes of radiation regions 
of all AGNs are larger than the so-called last stable radius 3R g , if we consider errors of 
calculated sizes. 



3.4. Accretion Rates 

The accretion rate M can be derived from the relation, L x ~ gm bhM 

Fix 

M = 2L x (R x /R g )/c 2 . (10) 

where L x is the soft X-ray luminosity in 0.1-2.4 keV band, G is gravitational constant, 
and R x is the size of the X-ray radiation region. Using R x /R g and L x , accretion rates for 37 
AGNs are obtained from Eq. 10, which are listed in column 6 of Table 1. The M distribution 
is < log(M) >= —1.50 ± 0.23 with a standard deviation of 1.42. In Figure 6 we show the 
accretion rates (M) versus black hole masses (Mbh)- The accretion rates of most of AGNs 
are 100-0.01 M EM , taking M Edd = 0.2^-iMe/yr) . 



We plot L x versus M in Figure 7. L x = (Mc 2 , where ( is the accretion efficiency The 
solid line in Figure 7 is L x = O.lMc 2 and the dash line is L x = O.OOlMc 2 . Most AGNs 
are in the range of 0.1~0.001 of (. Combining Eq. 10, we can drive the accretion efficiency 
( = 1/(2R X /R g ). From Figure 5 we find the values of R x /R g in most AGNs are in the range 
of 5~500, which is consistent with the value of ( (0.1~0.001). 
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3.5. The Compactness Parameter 

If the power law spectra observed in AGNs in the X-ray extend to much higher energies, 
then pairs may be expected to be produced if the optical depth of the 7-ray photons for 
pair production exceeds unity. The condition is usually stated in terms of the compactness 
parameter (Done & Fabian 1989), 

m e c 6 R x 

Where L x is the luminosity produced in the region of R x , m e is the mass of electron and cft is 
Thomson cross section. The values of the compactness parameters are listed in column 9 of 
Table 1. If I > 20n, the source has unity optical depth in X-ray for 7-ray absorption at about 
one Mev and a significant fraction of the source luminosity can then pass through electron- 
positron pairs. The distribution of the compactness parameter is < log(l) >= 0.90 ± 0.20 
with a standard deviation of 1.23. We plot the compactness parameter versus R x / Rg in 
Figure 8. A simple least square linear regression (Press et al. 1992) gives log(l) = (1.70 ± 
0.29) + (-0.85 ± 0.25)log(R x /R g ) (excluding PHL 1092) (R=-0.50, P=0.00168 ). Those 
AGNs with smaller R x /R g appear to have larger compactness parameter. From Figure 8 
we can find that several AGNs have / > 207r, while several AGNs have / < 0.1, which is 
consistent with the results of Done & Fabian (1989). 



4. THE SIZES OF THE BLR 

Broad emission lines are one of the dominant features of spectra of AGNs. Broad Line 
Regions (BLRs) play a particularly important role in our understanding of AGNs by virtue 
of their proximity to the central source. With reverberation mapping techniques, the sizes 
of the BLRs can be obtained through the study of correlated variations of the lines and 
continuum fluxes (Peterson 1993). The BLR sizes for 17 Seyfert 1 galaxies (Wandel et al. 
1999) and for 17 PG quasars (Kaspi et al. 2000) have been recently obtained. In Table 2, we 
list 12 AGNs with available BLRs sizes and the calculated accretion rates. Here we adopted 
3 light days as the error of the BLRs size of NGC3516, which is not given in Ho (1998). 
Errors of BLRs sizes of other 11 AGNs are from Kaspi et al. (2000). 

We assume that the gravitational instability of the standard thin disk at large radius 
leads to the formation of BLRs (Bian & Zhao 2002). The criterion is Q = < 1 (Golreich 
& Lynden-Bell 1965), where Q is the Keplerian angular velocity at R away from the central 
black hole (GM/R 3 ) 1/2 , M is the black hole mass, and p is the local mass density. We adopted 
the solutions of p for the standard thin disk from Shakura & Sunyaev (1973) , p = 3.1 x 
10~ 5 a~ 7 / 10 M^ 20 M^ 8 i?^ 4 15 ^ 8 f 11 / 5 g cm -3 , where a is the parameter of the standard a disk, 
/ = (1 - (R/R s ) 1/2 ) 1/4 , M 26 = M/(10 26 g s- 1 ), M 8 = M/(1O 8 M ), and R u = R/(10 14 cm). 
f = 1 and a = 1 are adopted in AGNs. We can obtain the sizes of the BLRs if we know the 
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central black hole masses and accretion rates, 

Ru = 880a 28 / 45 g- 8 / 9 M 2 - 22/45 M 8 1/3 . (12) 

The BLRs sizes of 12 AGNs are obtained from Eq. (12), which are listed in Table 2. 
The uncertainties of our calculated BLRs sizes are estimated by considering 1/4 < Q < 4 ( 
Golreich & Lynden-Bell 1965 ), errors of the black hole masses and the accretion rates. 

Figure 9 shows that the calculated BLRs sizes derived from Eq. 12 are consistent 
with the sizes from the reverberation mapping method. We have a x 2 (Press et al. 1992) 
test of logRsLR to show whether our calculated BLRs sizes are consistent with that from 
reverberation mapping method . x 2 an d probability are 0.933 and 99.996% for 12 AGNs. 



5. SUMMARY 

We use the exponential time-scale to characterize the soft X-ray variabilities of AGNs 
with obvious exponentially increasing or decreasing patterns in the light curves and obtain 
time-scales for 65 AGNs. The exponential time-scale is suitable for analyzing X-ray vari- 
ability, especially for the data with larger time gaps such as ROSAT data. The exponential 
time-scales can provide good upper limits for the central black hole masses and the sizes of 
the X-ray radiation region. From the relationship of L x and At e we can reject bremsstrahlung 
as the primary radiation source for most AGNs here. Many AGNs have very efficient engines. 
There is no explicitly grouped distribution in M — M diagram for different kinds of AGNs. 
This means that, under the same soft X-ray luminosity, the central mass or radiation region, 
compared with Schwarzschild radius, varies from one AGN to another. 

Using the recent mass estimates of central black holes and the exponential time-scales of 
AGNs from ROSAT/PSPC pointing data, the sizes of X-ray radiation region, the accretion 
rate and the compactness parameter of 37 AGNs are obtained. We use the calculated 
accretion rates to calculate the radius of gravitational instability for 12 AGNs among 37 
AGNs, which is consistent with the BLRs sizes from the reverberation mapping method. 
This provides further evidence for the gravitational instability leading to the formation of 
BLRs. 

We thank Keliang Huang for useful discussions, and the anonymous referee for the 
valuable comments. We thank Helmut Abt for reading our manuscript. We thank the 
Chinese Natural Science Foundation for financial support under contract 10273007. 
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Fig. 1. — Variability of count rate and exponential time-scale fitting of NGC 5548 in three 
energy band: 0.1-2.4 (total band), 0.1-0.4 (A band), 0.5-2.0 (B band) keV. The time zero is 
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Fig. 2. — The black hole masses versus the mass limits derived from the exponential time- 
scales. The solid line shows 1:1. 
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Fig. 3. — The soft X-ray luminosity versus the exponential time-scales. Solid line: Eddington 
limit; dash line: efficiency limit (r/=0.1) ; dot line: compact parameter limit (1=60); dash 
dot line: efficiency limit (r]—0. 007) ; dash dot dot line: bremsstrahlung limit, respectively, 
assuming L boi = 10L X . 
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Fig. 4. — The soft X-ray luminosity versus the black hole masses. The solid line shows 
L x = 0.1L Edd . and dash line is L x = 10~ 4 L EcW . 
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Fig. 5. — log(R x / Rg) versus log(MBu)- The dash line is R x /R g =l and the dot line is 
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Fig. 6.— log(M) versus log(M BH ). M = eM EDD . e is 1, 10, 0.1, 100, 0.01 respectively for 
the solid line, dash line, dot line, dash dot line and dash dot dot line. 
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Fig. 7. — log(L x ) versus log(M). The solid line and dash line show L x = (Mc 2 , £ = 0.1, 
C = 0.001, respectively. 
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Fig. 8. — Log of the compactness parameter I versus log(R x /R g ). The dash line is I = 20n. 
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Fig. 9. — Log(R ca i) versus log(RBLR)- The solid line shows 1:1. 
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Table 1: The properities of 65 AGNs in Right Ascension order. Col.l: name, Col. 2: type, 
Col. 3: the exponential time-scale, Col. 4: 0.1-2.4 keV Luminosity , Col.5:log of the BH mass 
in M , Col. 6: log of calculating accretion rates in M & /yr , Col. 7 log of the ratio of the 
size of X-ray emission region and the Schwarzschild radius , Col. 8: log of the compactness 
parameter. *: the exponential time-scale from the rising variation. Estimation of central 
black hole masses: a: Nelson 2000; b: Nelson & Whittle 1995; Falomo, et al. 2002; c: Wang 
& Lu 2001; d: Vestergaard 2002; e: Tremaine et al. 2002; f: Czerny et al. 2001; g: Mathur 
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Table 2: The BLRs size of 12 AGNs. Col.l: name, Col. 2: BLRs sizes from reverberation 
mapping method (Kaspi et al. 2000) , Col. 3: central black hole masses (Kaspi et al. 2000) , 
Col. 4: calculating accretion rate, Col. 5: calculating BLRs size. 



